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ABSTRACT 

We explore the sensitivity of the nucleosynthesis of intermediate mass elements (28 < 
A $80) in supernovae derived from massive stars to the nuclear reaction rates employed 
in the model. Two standard sources of reaction rate data (Woosley et al. 1978; and 
Thielemann et al. 1987) are employed in pairs of calculations that are otherwise identical. 
Both include as a common backbone the experimental reactions rates of Caughlan & 
Fowler (1988). Two stellar models are calculated for each of two masses: 15 and 25 
M . Each star is evolved from core hydrogen burning to a presupernova state carrying 
an appropriately large reaction network, then exploded using a piston near the edge of 
the iron core as described by Woosley & Weaver (1995). The final stellar yields from 
the models calculated with the two rate sets are compared and found to differ in most 
cases by less than a factor of two over the entire range of nuclei studied. Reasons for the 
major discrepancies are discussed in detail along with the physics underlying the two 
reaction rate sets employed. The nucleosynthesis results are relatively robust and less 
sensitive than might be expected to uncertainties in nuclear reaction rates, though they 
are sensitive to the stellar model employed. 
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1. INTRODUCTION 

Two major efforts to model nucleosynthesis in mas- 
sive stars have operated in parallel for a number of 
years (e.g., Woosley & Weaver 1995; Thielemann, 
Nomoto, & Hashimoto 1996). These groups have dif- 
fered in the stellar physics they employed - convection 
theory, treatment of stellar energy generation, start- 
ing point (the main sequence for the former, helium 
cores for the latter), and how the supernova explo- 
sion was simulated. They also differed, especially for 
elements heavier than silicon, in their choice of nu- 
clear reaction rates. Nomoto and his collaborators 
have always used the Hauser-Feshbach calculations of 
Thielemann et al. (1987; henceforth TAT). So far, 
Woosley & Weaver have used those of Woosley et al. 
(1978; henceforth WFHZ). 

Considering the many differences between the mod- 
els, it is comforting (perhaps even surprising) that 
the predicted nucleosynthetic yields from these two 
groups agree as well as they do, but detailed com- 
parison does reveal discrepancies. Here we shall seg- 
regate those differences that exist because of the di- 
chotomy of stellar models from those reflecting purely 
the choice of nuclear physics. We shall also discuss the 
underlying assumptions in each reaction rate set and 
why the two differ. 

We begin by describing (§2) the two approaches 
to stellar evolution and, briefly, the two rate sets 
(pending a more thorough review in §4). Section 3 
presents new yields for some hybrid models - stars 
evolved using the KEPLER code of Weaver, Zimmer- 
man, & Woosley (1978) and exploded as described in 
Woosley & Weaver (1995), but using rates from both 
groups (TAT & WFHZ). In addition to helping to 
understand why calculations of the two groups differ, 
the use of two independent sets of reaction rates in 
identical stellar models helps determine the (nuclear 
physics portion of the) error bar one should assign 
to nucleosynthesis studies of this sort. We find, in 
general that it is relatively small. 

In §4, we examine in greater detail the nuclear 
physics that produced the two reaction rate sets, 
specifically the input data and physical assumptions 
employed in two globally parameterized Hauser-Feshbach 
calculations by Woosley et al. (1978) and Thielemann 
et al. (1987). Previously unpublished comparisons to 
experimental rates are given along with mutual com- 
parisons of the two sets. 

Section 5 then discusses the differences in stellar 



yields in light of the results of §4 as well as the known 
differences in the stellar models. In §6, we summarize 
and conclude. 

2. Previous Models and the Rates They Em- 
ploy 

2.1. A Tale of Two Stellar Models 

The nucleosynthesis calculations provided to the 
community by Thielemann, Nomoto & Hashimoto 
(1996, hereafter TNH) are based on the stellar evo- 
lution calculations of Nomoto & Hashimoto (1988) 
for massive stars. These calculations were extended 
to a larger grid of progenitor masses by Hashimoto 
et al. (1993). Both works followed the evolution of 
helium cores of fixed mass (omitting hydrogen burn- 
ing and complications of mass loss and red giant for- 
mation) and both used the Schwarzschild criterion 
for convective instability. The nuclear energy gen- 
eration and nucleosynthesis within a structural iter- 
ative loop (stellar evolution time-step) was evaluated 
using a reduced reaction network adopted to the ap- 
propriate burning stages (for details see Thielemann, 
Hashimoto, & Nomoto, 1990; Nomoto & Hashimoto 
1988; Hashimoto, Hanawa, & Sugimoto 1983; and 
Hashimoto, 1995). A 30 isotope network was used 
from He to oxygen burning that included essentially 
only stable nuclei up to Ca and did not include 
neutron-rich nuclei. Silicon burning, a much more 
complicated stage of stellar burning in which equili- 
brated strong flows must be followed and weak inter- 
actions on trace isotopes are important, was treated 
in a different way. During the onset of silicon burning 
a 148 isotope network was employed, extending to Ca 
and including also unstable and neutron-rich species. 
After partial burning of Si, when quasi-equilibrium in 
the Si-group was attained, the calculation switched to 
using tables of the relevant quantities evaluated previ- 
ously with a 299 isotope network given in Thielemann, 
Hashimoto, & Nomoto (1990). The tables gave for 
each instant of silicon burning at density, p, temper- 
ature, T, and electron mole number, Y e , the appropri- 
ate energy generation rate, silicon depletion rate, and 
time derivative of Y e . Technically this was done by 
interpolating quantities from prior calculations for a 
grid of densities, temperatures, Y e , and Si mass frac- 
tion. The individual abundances of the isotopes could 
be obtained, if necessary, from quasi-equilibrium at 
the same p, T, Xgi, and Y e . 

Such a treatment led to an approximately correct 
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energy generation rate during all stages of stellar evo- 
lution and also kept track of the major abundances. 
Limitations in this network, however, which e.g. did 
not contain all production or depletion reactions of 
26 Al, meant that it could not give reliable answers 
for the abundances of such isotopes. Similar limi- 
tations on the neutron-rich side of the (stellar evo- 
lution) network did not permit the s-process to be 
calculated. "Post-processing calculations" by Prant- 
zos, Hashimoto & Nomoto (1990) gave an indication 
of what s-processing might be expected during He- 
burning, but this study centered on the major s- 
process nuclei beyond Fe (and not lighter elements 
like Si, P, S, CI, Ar, K, Ca, and Ti as discussed in 
Kappeler et al. 1994) 

This treatment differed in many ways from Wooslcy 
& Weaver (1995, hereafter WW95). These authors 
also made use of a smaller 19 isotope network to 
calculate energy generation during hydrogen, helium, 
carbon, neon, and part of oxygen burning. Though 
the matrix contained only 19 species, assumed steady 
state links for about a dozen other nuclei made it 
the equivalent of a larger network of about 30 iso- 
topes (Weaver, Zimmerman, & Woosley 1978). After 
the oxygen abundance declined, during central oxy- 
gen burning, to below 4%, the transition was made 
to a larger quasi-equilibrium network of 137 isotopes. 
This network was converged during each iterative cy- 
cle of the stellar code and accounted for most of the 
computational expense of the entire evolutionary cal- 
culation. In addition, Weaver et al. (and subse- 
quently WW95) used a different criterion for con- 
vective instability, the Ledoux criterion plus semi- 
convection in the regions that were Schwarzschild un- 
stable but not Ledoux unstable. This contributed sig- 
nificantly to the variable iron core masses obtained 
by the two groups for a given helium core mass (in 
the Ledoux criterion a Y e barrier inhibits convection, 
in the Schwarzschild criterion it does not). Finally, 
WW95 treated their large network nucleosynthesis in 
a different way. A reaction network of 200 isotopes 
(large by the standards of the time) was carried in a 
"co-processing" mode in each zone throughout the en- 
tire stellar evolution. Thus the converged T and p of 
each stellar model was used to update this network at 
each time step, and most importantly, convective mix- 
ing was carried out in the co-processed abundance set 
wherever appropriate. The s-process up to A=65 was 
also calculated. This sort of detail was missing in the 
presupernova models of Nomoto & Hashimoto (1988) 



and Hashimoto et al. (1993), and thus restricted the 
information available from the presupernova evolu- 
tion when it came time to calculate the explosion. 

During the explosive part of the modeling, there 
were differences as well, though not so dramatic. The 
size of the nuclear networks for this phase have been 
discussed in WW95 and in Thielemann, Hashimoto & 
Nomoto (1990). Because of persistent uncertainties in 
modeling the explosion of Type II supernovae as well 
as the practical need to explore a wide range of stellar 
models, both groups simulated explosions artificially. 
WW95 made use of an inward (during in-fall, before 
the explosion) and outward (after the induced explo- 
sion) moving piston with a velocity which caused a 
final kinetic energy of the ejecta of 1.2 — 2 x 10 51 erg. 
The details of this method were first described in 
Woosley & Weaver (1986) and are discussed in de- 
tail in WW95. Thielemann, Hashimoto & Nomoto 
(1990), Hashimoto et al. (1993), TNH, and Nomoto 
et al. (1997) made use of a "thermal bomb", i.e. in- 
troducing thermal energy in the form of a tempera- 
ture enhancement inside the Fe-core of the progenitor 
star. The thermal energy was also chosen in such a 
way to attain a total kinetic energy of the ejecta of 
10 51 erg. While WW95 obtained a mass cut between 
ejecta and neutron star (which determined the total 
amount of 56 Ni ejected) from their choice of piston po- 
sition and energy (though this mass cut often lay well 
outside their piston), TNH adjusted their mass cut 
based upon existing observations which relate super- 
nova progenitor masses to the 56 Ni ejecta. Both treat- 
ments have their advantages and disadvantages (see 
Aufderheide, Baron, & Thielemann 1991), neither is 
an adequate substitute for a full understanding of the 
supernova mechanism. They are parameterizations, 
but the parameters do reflect real physics, motion in- 
wards as well as outwards for the piston model, high 
entropy for the innermost ejecta (reflecting the blow- 
ing of a bubble by neutrino energy deposition) in the 
bomb model. 

2.2. The Two Reaction Rate Libraries 

In addition to the differences in stellar physics em- 
ployed by both groups, each employed an indepen- 
dent library of nuclear reaction rates. These two data 
sets are a compilation of both experimental and theo- 
retical reaction rate information. The rate sensitivity 
problem separates nicely into two mass regions, above 
and below silicon. 

Below silicon, most of the important rates deter- 
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mining energy generation and nucleosynthesis up to 
core carbon burning have been measured directly in 
the laboratory [Caughlan & Fowler (1988, hereafter 
CF88); Caughlan et al. (1985); Harris et al. (1983); 
and Fowler, Caughlan, & Zimmerman (1975)]. Fits 
to these reaction rates have been adopted for some 
time by both groups and form a common backbone of 
the two independent reaction rate libraries with one 
important exception, the rate adopted for the critical 
reaction 12 C(a,7) 16 0. WW95 used a fit to this reac- 
tion rate recommended by CF88 but multiplied at all 
temperatures by a constant factor of 1.7. This corre- 
sponds to an S- factor (at 300 keV) of 170 keV barns, a 
value which has been determined as optimal for pro- 
ducing the solar abundance set (Weaver & Woosley 
1993; Timmes, Woosley k Weaver 1995). Nomoto & 
Hashimoto (1988) chose the rate for 12 C(a,7) pre- 
scribed by Caughlan et al. (1985), which corresponds 
at T g = 0.3 (i.e. core helium burning) to CF88 * 2.35. 
Both of these values are compatible with the still ex- 
isting uncertainties of this rate (Buchmann 1997). An 
examination of the impact on nucleosynthesis in the 
stellar models due to the particular prescription for 
this critical rate is deferred to §3. 

In the heavier mass range, reaction rates, except 
for the relatively well known (n,7) cross sections (Bao 
& Kappeler 1987; Wisshak et al. 1997), are chiefly 
a product of Hauser-Feshbach theory. For the the- 
oretical rates, there have basically been two choices 
available - we present calculations using both, which 
will be loosely referred to as "TAT rates" (Thiele- 
mann, Arnould, & Truran 1987) and the 'WFHZ 
rates" (Holmes et al. 1976; Woosley et al. 1978). 
The corresponding Hauser-Feshbach codes that cal- 
culated the nuclear cross sections from which these 
rate sets were constructed are loosely referred to as 
"SMOKER" (TAT) and "CRSEC" (WFHZ, Holmes 
1976). 

The neutron capture reaction rates in both reac- 
tion rate libraries incorporate the experimental cross 
section recommendations of Bao & Kappeler (1987, 
for a more recent compilation see Beer, Voss, & Win- 
ters, 1992). Above silicon, both groups have sup- 
plemented their Hauser-Feshbach reaction rates with 
results from experiment. The TAT set included re- 
actions on neutron- and proton-rich unstable tar- 
gets from Malaney & Fowler (1988, 1989), Wiescher 
et al. (1986, 1987, 1989), and Van Wormer et al. 
(1994). Beyond this the TAT reaction rate library 
depended strictly on rates derived from the Hauser- 



Feshbach calculations of Thielemann, Arnould, and 
Truran (1987). 

For charged-particle induced reactions on interme- 
diate mass nuclei, the experimental reaction rates im- 
plemented in the WFHZ reaction library are mostly 
drawn from D.G. Sargood and his collaborators, Mitchell, 
Tingwell, Tims, Hansper, and Scott. In most in- 
stances, this group has generated fits to reaction rates 
derived from their exhaustive cross section measure- 
ments, these have been adopted in the WFHZ set 
wherever given. When such fits were not available 
(Sargood 1982), the fits to the reaction rates gen- 
erated by the Hauser-Feshbach studies of Woosley 
et al. (1978) were normalized to agree with ex- 
periment at temperatures appropriate to the nuclear 
burning stage where the reaction rate was most im- 
portant. For rates proceeding on targets heavier than 
Si, this was usually Tg — 3.0 (i.e. oxygen and sili- 
con burning). Comparisons between experiment and 
the Hauser-Feshbach results of both Thielemann et 
al. (1987) and Woosley et al. (1975, 1978) are pre- 
sented in §4. A complete bibliography of all the 
rates employed and the authors contributing to the 
WFHZ reaction rate compilation used in this study, 
including an extensive tabulation of these rates for 
0.03 < Tg < 10, is given in Hoffman & Woosley 
(1992). This data, plus the fit constants for the 
TAT rates, are also available on the world wide web 
( bttp: / /ic. lbl.gov/astro /astroratc.html ) . 

For weak reactions, both groups always included as 
lower bounds to the weak decay rates the ground-state 
rates taken from the measured laboratory values given 
in the Nuclear Wallet Cards (1990). Also included in 
all calculations were electron captures on protons and 
positron captures on neutrons, which proved impor- 
tant. For temperatures above 5 x 10 s K and den- 
sities greater than 10 5 g cm -3 , all models used the 
weak interaction rates of Fuller, Fowler, & Newman 
(1980, 1982a,b, 1985) for all nuclei lighter than mass 
60. Above mass 60, no reliable tabulation of the tem- 
perature and density dependence of the /3-decay rates 
are available, the ground state rates were used. 

Finally, neutrino induced reactions on nuclei were 
included in the manner described in WW95. These 
are based on the work of Wick Haxton as described in 
Woosley et al. (1990). For select neutrino energies of 
4.0, 6.0, & 8.0 MeV, tables of cross sections and par- 
ticle evaporation branching ratios for charged-current 
and neutral-current inelastic neutrino scattering reac- 
tions on intermediate mass nuclei are given by Hoff- 
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man & Woosley (1992). Neutrino reaction rates have 
not been previously incorporated in either TNH or 
Nomoto et al. (1997), but are included in all models 
presented in this work. 

3. New Hybrid Results 

To explore the dependence of the nucleosynthesis 
strictly upon the choice of reaction rates above A=28, 
a set of new hybrid models was generated that used 
the Woosley- Weaver choice of stellar physics, but var- 
ied the reaction rate library. Two stars having main 
sequence masses of 15 and 25 M were evolved from 
core hydrogen burning to presupernova collapse and 
subsequent explosion. Each calculation was repeated 
twice, once with each rate set. Since the network em- 
ployed for nucleosynthesis is decoupled from that used 
for energy generation and Y e changes in the star, it 
was feasible to split these operations. 

This also means that the energy generation net- 
work stayed the same in both cases (including the 
nuclear input), which guaranteed that the tempera- 
ture and density conditions between two models of 
the same mass but using different reaction rate sets 
stayed exactly the same, thus only probing the rate 
uncertainties for the (co-processing) nucleosynthesis. 
However, this did lead to a minor inconsistency be- 
tween the rates used in the small energy generation 
network used to converge the stellar model and their 
counterparts in the large TAT reaction rate library 
utilized for nucleosynthesis. For all reactions on tar- 
gets lighter than silicon, the rates in both sets and in 
the energy generation network are drawn from CF88, 
including the rate for 12 C (a, 7) 16 (multiplied by 1.7) 
and the heavy ion-rates ( 16 0+ 16 0, etc.). For heavier 
species contained within the energy generation net- 
work, the Hauser-Feshbach rates of Woosley et al. 
(1978) were used. This will have little effect on the re- 
sults of the converged stellar model, in that the struc- 
ture of the star is mostly determined prior to oxygen 
burning, where the CF88 rates will play the major 
role. 

3.1. Stellar Yields for the WFHZ and TAT 
Reaction Rate Libraries 

Figures |l| and || compare the detailed stellar nu- 
cleosynthesis from four sets of calculations. Plotted 
for each stable isotope (after decay) is the ratio of 
the stellar yield (in solar masses) from the models us- 
ing the "TAT rates" divided by the equivalent num- 



ber obtained using "WFHZ rates" . The stellar yields 
from the calculations using the WFHZ rate set should 
have been, and in general are equivalent to those from 
models S15A and S25A of WW95. The networks em- 
ployed in each survey extended to 90 Kr, (see Table 
|l|), although the range of nuclei shown in the figures 
only extends to 76 Ge. Isotopes of the same element 
are connected by solid lines. The most abundant iso- 
tope of a given element is plotted as an asterisk. Nu- 
clei made chiefly as radioactive progenitors are sur- 
rounded by diamonds, squares, or circles, depending 
on whether the given nucleus was made by the same 
or different radioactive species as a result of the dif- 
ferent rate sets used in the model studied (if a nu- 
cleus was chiefly produced as the same radioactive 
progenitor by both reaction rate sets, the choice was 
a diamond) . The dotted lines indicate a factor of two 
difference in the ratio of the stellar yield between the 
two models. Overall the comparison is quite good, 
much better than a factor of two in almost all cases. 
As we shall see (§4), the actual difference in reac- 
tion rates is frequently larger than a factor of two, 
but there is compensation, in that the major flows 
proceed through channels within the valley of beta 
stability that are in good agreement between the two 
rate compilations. 

The nuclei which show greater than 20% deviations 
in the 15 M and 25 M Q stars are: 33 S, 40 Ar, 40 K, 

44,46 Cai 45 SC) 50 Vj an d 70 Zn The causes for the£ , e 

discrepancies are analyzed is §5. 

3.2. Comparison of Stellar Yields to TNH 

When considering a direct comparison of stellar 



yields between the two groups (Woosley & Weaver 



1995; Thielemann, Nomoto, fc Hashimoto 1996| ), apart 
from the many differences in the separate treatments 
of stellar evolution (§2), an important distinction is 
the assumed rate for 12 C(a,7) 16 0. The rate for this 
reaction remains uncertain and yet plays a major role 
in determining both the final structure of the presu- 
pernova star and its ultimate nucleosynthesis (Weaver 
& Woosley 1993). TNH used the rate recommended 
by Caughlan et al. (1985), while WW95 adopted the 
rate recommended by Caughlan & Fowler (1998), but 
multiplied at all temperatures by a constant factor of 
1.7. In order to make a meaningful comparison be- 
tween these two models, as well as to explore the effect 
of this critical rate, the 15 M model of the previous 
section was recalculated using both reaction rate sets, 
but assuming a constant multiplicative factor for the 
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CF88 12 C(a,7) rate of 2.35, which agrees with the 
Caughlan et al. (1985) rate at Tg = 0.30 during core 
helium burning. This facilitates comparison between 
stellar models that otherwise would have been very 
different. 

The results are presented in Table ||. For a solar 
metallicity 15 M star, three sets of stellar yields (in 
M Q ) are given, those calculated with the WFHZ rate 
set, those with the TAT rate set, and finally the yields 
from the 15 M Q star of TNH. This list covers the sta- 
ble isotopes included in the TNH hydrostatic network 
up to oxygen, and thereafter all stable isotopes (after 
decay) up to 70 Ge. Also given are two ratios of the 
stellar yields, which compare yields using the same 
stellar model (WFHZ/TAT), and a common set of nu- 
clear reaction rates (TNH/TAT) in explosive burning. 
Blanks appear in place of entries that either were not 
calculated by TNH or for which the comparison was 
less than 0.001. The last column indicates the abbre- 
viation for the dominant nucleosynthetic process(es) 
responsible for production of the quoted stellar yield 
(or ratio of yields) for a given isotope (see Table 19 
WW95, also Woosley & Weaver 1999). 

Differences of more than 10% are uncommon be- 
tween calculations that use a common stellar model 
(WFHZ/TAT), but larger differences are apparent if 
the model varies (TNH/TAT). For the alpha-isotopes 
up to 28 Si the ratios are within a factor of three (typ- 
ically better) and are most likely due to the stel- 
lar model treatment (with the different treatments 
of convection playing the dominant role). Some 
agree remarkably well ( 20 Ne, 26 Mg, 31 P, 44 Ca, 48 Ti, 
60 ' 61 ' 62 Ni), while others are very discrepant ( 19 F, 
46 > 48 Ca, 50 Ti, 54 Cr, 58 Fe, and 64 Ni). These differ- 
ences reflect either the physics inherent to each stel- 
lar model or the limited reaction network used during 
stellar evolution phases in TNH (discussed in §2). 

14 N and 23 Na, for example, are products of hydrogen- 
burning, which is not calculated in TNH, they started 
with helium cores. The smaller differences within a 
factor of a few for nuclei close to the stable N = 
Z line up to calcium are most likely due the stel- 
lar model treatment. 19 F (and to a lesser degree 
26 Al) are discrepant because TNH did not imple- 
ment the ^-process (Woosley et al 1990). The lim- 
ited network used by TNH during hydrostatic stel- 
lar evolution likely accounts for the low production of 
many neutron-rich stable nuclei processed by neutron- 
capture during the s-process ( 36 S, 40 Ar, 46 Ca, 58 Fe). 
Yet other discrepancies are due to their limited pro- 



duction in massive stars. Ca, 50 Ti, and Cr, for 
example, are produced in nuclear statistical equilib- 
rium in SN Type la (Woosley & Eastman 1995). The 
yields quoted for these species in the Kepler based 
models resulted from s-processing. 

For nuclei predominantly made during explosive 
oxygen and silicon burning, the agreement is often 
good (x4). An exception are most of the Cl-Sc iso- 
topes, which typically differ by factors of 10 or more. 
These are susceptible to convective processing in the 
oxygen shell shortly before collapse. The s-process 
can also contribute to species in this mass range. 

Table |^ gives the important gamma-line radioac- 
tivities made in each model; 22 Na, 26 Al, 44 Ti, 56 Ni, 
60 Fe and 60 Co. (The quoted stellar yields have been 
included in the appropriate stable isobars in Table ||). 
Production of 44 Ti and 56 Ni agree within a factor of 
two, indicating that the different treatments affect- 
ing the mass-cut (moving piston vs. thermal bomb) 
achieve (in this case) consistent agreement between 
these two important supernova diagnostics. The oth- 
ers disagree for reasons stated above. 60 Fe and 60 Co 
differ because they were made via s-processing dur- 
ing neon-burning and in explosive processing in the 
Kepler based models. The 22 Na agreement is good, 
but the yield is very small. We think this isotope is 
made predominantly in novae. Convective process- 
ing in the neon-oxygen shell prior to collapse (in the 
Kepler models) produced the larger 26 Al yields. 

We emphasize that the choice of reaction rates 
above silicon have less of an effect on nucleosynthesis 
in massive stars than the stellar environment. Im- 
provements to both are clearly needed and welcome. 

4. A Comparison between Reaction Rate Sets 

4.1. Theoretical Nuclear Reaction Cross Sec- 
tions 

We now move to a discussion of the nuclear physics 
employed in the two studies. The reaction rate sets 
employed by both contain, as backbones, a number 
of experimentally determined values, the most impor- 
tant being the charged-particle rates of Caughlan & 
Fowler (1988) and the neutron capture rates of Bao 
& Kappeler (1987). The two sets differ however in 
the values adopted for the large number of reactions 
that either cannot, or have not been measured in 
the laboratory. A traditional theoretical approach in 
this mass range is the statistical or Hauser-Feshbach 
model. This model is valid only for high level densi- 
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ties in the compound nucleus, a condition that is gen- 
erally (but by no means universally) satisfied during 
the advanced evolution of massive stars because of the 
heavy fuels and high temperatures involved. The two 
compilations used respectively the Hauser-Feshbach 
calculations of Holmes et al. (1976) and Woosley et 
al. (1978) (WFHZ), or of Thielemann et al. (1987) 
(TAT). It is therefore advisable to have a short look 
at the different approach of these two nuclear models 
in order to understand the different results. 

A high level density in the compound nucleus al- 
lows one to use energy averaged transmission coeffi- 
cients T, which describe absorption via an imaginary 
part in the (optical) nucleon-nucleus potential (for de- 
tails see Mahaux and Weidenmuller 1979). This leads 
to the well known expression 



(2Jf + l)(2J i + l) 



E 



(2 J + 1) 



J,7T 



T tot (E, J, 7r) 

x T?(E,J,*,E?,J?,k?) 
x T»{E,J,v,E v m ,r m ^ m ) (1) 

for the reaction i^fj, o)m L ' from the target state i M to 
the exited state m v of the final nucleus, with center 
of mass energy E i3 and reduced mass fXij . J denotes 
the spin, E the excitation energy, and 7r the parity of 
excited states. When these properties are used with- 
out subscripts they describe the compound nucleus, 
subscripts refer to the participating nuclei i and m 
and projectile and emitted particle j and o in the re- 
action o)m u , where the superscripts indicate the 
specific excited states in the target nucleus i and final 
nucleus m. Experiments measure ^ v ^ u (jiO]Eij), 
summed over all excited states of the final nucleus, 
with the target in the ground state. Target states p 
in an astrophysical plasma are thermally populated 
and the astrophysical cross section <j*(j,o) is given 

by 

0-* (j, °; En ) = E a i U o; E v ) 

v 

x YP J i +l)exp(-£f/fcT) 

x ^(2Jf + l)cxp(- J Ef/fcT)^(2) 
The summation over v replaces T„(E, J, w) in Eq.([j]) 



by the total transmission coefficient 
T (E,J,w) = J2T^(E,J,n,E^,J^,w^) 



v=0 
E-S„ 



J ^2 T (E,J,n, E m , J m , 7r m ) 



x p(E m ,J m ,Tr m )dE„ 



(3) 



Here S m ^ is the channel separation energy, and the 
summation over excited states above the highest ex- 
perimentally known state v m is changed to an inte- 
gration over the level density p. The summation over 
target states /i in Eq. (||) has to be generalized accord- 
ingly. 

In addition to the ingredients required for Eq.(Q), 
like the transmission coefficients for particles and pho- 
tons or the level densities, width fluctuation correc- 
tions W(j, o, J, 7r) have to be employed as well. They 
define the correlation factors with which all partial 
channels of incoming particle j and outgoing parti- 
cle o, passing through excited state (E, J, n), have to 
be multiplied. This is due to the fact that the de- 
cay of the state is not fully statistical, some memory 
of the process of formation is retained and influences 
the available decay choices. The major effect is elastic 
scattering, the incoming particle can be immediately 
re-emitted before the nucleus equilibrates. Once the 
particle is absorbed and not re-emitted in the very 
first (pre-compound) step, equilibration is very likely. 
This corresponds to enhancing the elastic channel by 
a factor Wj. Besides elastic scattering, the main ef- 
fect is a redistribution of widths into weak channels, 
having a major impact at channel opening energies 
when i.e. a (p/y) and a (p,n) channel compete and 
the weaker (p,7) channel is enhanced. Verbaatschot, 
Weidenmuller and Zirnbauer (1984) derived a general 
expression in closed form, which is however computa- 
tionally expensive to use. A fit to results from Monte 
Carlo calculations was given by Tepel et al. (1974). 
While the width fluctuation corrections of Tepel et al. 
(1974) are only an approximation to the correct treat- 
ment, Thomas, Zirnbauer & Langanke (1986) have 
shown that they are quite adequate. For a general 
discussion of approximation methods see Gadioli and 
Hodgson (1992) and Ezhov and Plukjo (1993). Of 
the presently available thermonuclear rate sets in the 
literature (and used for nucleosynthesis calculations) 
this effect is only included in the TAT rate set. 

The important ingredients of statistical model cal- 
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culations are the particle and 7-transmission coeffi- 
cients T and the level density of excited states p. 
Therefore, the reliability of such calculations is de- 
termined by the accuracy with which these compo- 
nents can be evaluated. In the following we discuss 
the methods utilized to estimate these quantities in 
the two sets of reaction rates. 

4.. 1.1. Particle Transmission Coefficients 

The transition from an excited state in the com- 
pound nucleus (E,J,n) to the state (E'f, Jf,7rf) in 
nucleus i via the emission of a particle j is given by 
a summation over all quantum mechanically allowed 
partial waves 



T?(E,J,*,E?,J?,i$) 



J i+ J i 



J+s 

l=\J—s\ S=\J^-Jj\ 



T JlB (Et 



(4) 

Here the angular momentum I and the channel spin 
s = Jj + couple to J = I + s. The individual 
transmission coefficients Ti are calculated by solving 
the Schrodinger equation with an optical potential 
for the particle-nucleus interaction. The studies of 
thermonuclear reaction rates by Holmes ct al. (1976) 
and Woosley et al. (1978) employed optical square 
well potentials and made use of the black nucleus 
approximation. Thielemann et al. (1987) used the 
optical potential for neutrons and protons given by 
Jeukenne, Lejeune, & Mahaux (1977), based on mi- 
croscopic calculations, with the local density approx- 
imation. It included the corrections of the imaginary 
part by Fantoni, Friman, & Pandharipande (1981) 
and Mahaux (1982). For a-particles they employed a 
phenomcnological Woods-Saxon potential derived by 
Mann (1978), based on extensive data by McFadden 
& Satchler (1966). 

The differences between these treatments can be 
seen in Figures ||, |J, and ||, which show the ratios of 
the TAT/WFHZ s-wave transmission coefficients for 
a set of stable nuclei with varying mass number A at 
a number of particle energies (the large jumps in the 
ratios over a narrow range in A are due to this choice 
of abscissa, with increasing mass one jumps from iso- 
tope to isotope as well as element to element). The 
ma in differences arc due to the black nucleus approx- 



irn qtion made by WFHZ. This is equivalent to a fully 



absorptive potential, once a particle has entered the 
potential well, and therefore does not permit reso- 
nance effects. If single particle s-wave states are lo- 



cated close to the top of the potential well, they can be 
experienced as resonances by incoming (s-wave) par- 
ticles. This is shown for neutrons in Figure |3|. The 
maxima occur for compound nuclei with 2si/ 2 and 
3si/2 configurations at mass numbers Am50 and 150. 

Increasing particle energies leads to penetration of 
the potential well at energies further removed from 
the resonance energy for the same target nuclei. Thus, 
the difference to the black nucleus approach, for the 
nuclei experiencing the resonance behavior at low en- 
ergies, vanishes and the ratio approaches unity. We 
see, however, that for thermal neutron energies of 
0.1 MeV (T w 10 9 K), the difference is still quite re- 
markable (dotted line, Figure^). Lighter nuclei will 
have the same (unpopulated) 2si/ 2 and 3si/ 2 configu- 
rations at higher excitation energies in the continuum, 
where they lead to very broad resonances and there- 
fore the ratio also approaches unity. This changes 
for p-wave neutrons, interacting with lpi/2.3/2 and 
2pi/2,3/2 configurations, because the centrifugal bar- 
rier can still cause sharper resonances in the contin- 
uum, and with increasing energy a shift of the peaks 
(which however decline in size) to lower mass numbers 
is observed (not shown here). 

Proton s-wave states below the Coulomb barrier, 
can, similar to the centrifugal barrier behavior for p- 
wave neutrons discussed above, exist as resonances 
for unbound particle energies (similarily the peaks 
shift with bombarding energy). Therefore, strong res- 
onance absorption still occurs for proton energies of 
2-5 MeV (see Fig. ||). In astrophysical applications 
these are energies in the Gamow window, which cor- 
respond to temperatures of 2-5 xl0 9 K. 

An additional effect, which is only pronounced for 
a-particles, can be seen in Fig. §. The ratio of the 
transmission functions is systematically rising with 
mass number. For lower a-energies the rise starts 
at lower mass number. This behavior can be under- 
stood in terms of absorption in the Coulomb barrier 
( [Michaud, Scherk, fc Vogt 1970p. At sub-Coulomb 



energies the particle wave function is increasing ex- 
ponentially as a function of radius within the bar- 
rier. This has the effect that important contributions 
to the transmission functions occur well beyond the 
normal nuclear radius at energies very far below the 
Coulomb b arrier (cf. Fig. 11 of Michaud, Scherk, fc 
Vogt 1970). These contributions are not considered in 



the treatment of square well potentials and the black 
nucleus approximation, which assumes that absorp- 
tion takes place inside the nuclear radius. Therefore, 



8 



the transmission is systematically underestimated for 
sub- Coulomb energies. 

It should be noted, however, that a correction fac- 
tor / is introduced within the framework of square 
well potentials, which corrects for the unphysical ad- 
ditional reflection at the edge of the potential well. It 
is taken to be independent of energy which has been 
shown to be a good approximation for nucleonic pro- 
jectiles (Peaslee 1957, Vogt 1968). However, it should 
be energy-dependent for a-particles. It is an interest- 
ing fact that the two approximations made in neglect- 
ing absorption deep in the barrier and in employing 
an energy-independent correction factor nearly cancel 



for energies not too much below the barrier ( Michaud 
Scl jerk, fc Vogt 1970[ ). That is why the ratio of 
the TAT and WFHZ a-transmission functions devi- 
ate from unity only for A > 60, even for low energies. 
(Historically, this is also the reason why an energy- 
dependent correction factor for a-particles was rarely 
used because it did not improve the accuracy of the 
transmission functions.) 

4-1-2. "/-Transmission Coefficients 

The dominant 7-transitions (El and Ml) have to 
be included in the calculation of the total photon 
width. The smaller, and therefore less important, Ml 
transitions have usually been treated with the simple 
single particle approach (T cx E 3 , Blatt and Weis- 
skopf 1952), as discussed in Holmes et al. (1976) 
and used in both compilations presented here. The 
El transitions are usually calculated on the basis 
of the Lorentzian representation of the Giant Dipole 
Resonance (GDR) [see the experimental compilations 
which are provided e.g. by Berman (1975), Berman 
and Fultz (1975), Carlos et al. (1974), Berman et al. 
(1979), and Gurevich et al. (1981)]. The El transmis- 
sion coefficient for the transition emitting a photon of 
energy E 1 is a function of E 1 , the GDR energy Eq 
and width Tq. Holmes et al. (1976) and Woosley et 
al. (1978) used analytical fits as a function of A and 
Z which are quite accurate for Eq. The (hydrody- 
namic) droplet model approach for Eq by Myers et 
al. (1977) gives an excellent fit to the GDR energies 
and can also predict the split of the resonance for de- 
formed nuclei, when making use of the deformation. 

The width of the GDR is less understood even in 
modern Random Phase Approximation calculations 
(Catara et al. 1997), especially for the pronounced, 
observed shell effects with strongly reduced widths at 
magic numbers. Holmes et al. (1976) and Woosley et 



al. (1978) used an analytical fit to simulate this be- 
havior. Thielemann and Arnould (1983) proposed a 
phenomenological model for the GDR width which 
satisfactorily reproduces the experimental data for 
spherical and deformed nuclei and can be described 
as a superposition of a macroscopic width due to 
the viscosity of the nuclear fluid and a coupling to 
quadrupole surface vibrations of the nucleus (see also 
Cowan, Thielemann & Truran 1991). 

Utilizing the methods outlined above for the to- 
tal radiation width at the neutron separation energy 
for all nuclei with experimental information, Thiele- 
mann et al. (1987) found agreement generally within 
a factor of 1.5. This was an improvement over earlier 
attempts by Johnson (1977) and Hardy (1982). We 
expect from several tests that the differences between 
the WFHZ and TAT gamma transmission functions 
is not larger than a factor of 2. 

4.1.3. Level Densities 

Most statistical model calculations use the level 
density description of the back-shifted Fermi gas (Gilbert 
and Cameron 1965) 



1 



p{U,J,n) = -f(U,J)p(U) 
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which assumes that positive and negative parities 
are evenly distributed and that the spin dependence 
f(U, J) is determined by the spin cut-off parameter a. 
The level density of a nucleus is therefore dependent 
only on two parameters: the level density parameter 
a and the back-shift S. Gilbert and Cameron (1965) 
were the first to identify an empirical correlation with 
experimental shell corrections S(N, Z) 



- = c + ciS{Z,N), 
where S(N, Z) is negative near closed shells. 



(7) 
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Holmes et al. (1976) made an extensive tabula- 
tion of a/ A as a fit dependent on N and Z in different 
nuclear mass ranges, which can be used to predict 
the level density for any stable and unstable nucleus. 
Woosley et al. (1978) also made use of a large amount 
of experimental level density information, available 
for stable nuclei. This resulted in a set of experimen- 
tal level densities for such nuclei, which was utilized 
in their statistical model calculations. 

Thielemann et al. (1987) kept the same formalism 
as in Eq. (^) but performed an independent evaluation 
of the coefficients cq and c\. By dividing the nuclei 
into three classes (a: those within three units of magic 
nuclcon numbers, b: other spherical nuclei, c: de- 
formed nuclei), an improved agreement with respect 
to all previous approaches was obtained (maximum 
deviations of a factor of 3 to 5 at the neutron sepa- 
ration energy) . This treatment was very phenomeno- 
logical and unsatisfactory and still the weakest point 
in their procedures for the calculation of cross sec- 
tions. Further improvements of level density predic- 
tions have been recently attained (Rauscher, Thiele- 
mann & Kratz 1997), but are not yet included in the 
rate sets which are compared in the present paper. In 
general deviations in level densities of up to a factor 
of 2-3 (if not more) have to be expected for the two 
rate sets. 

Reverse rates and rates of thermally populated tar- 
gets include partition functions of the participating 
nuclei. These depend on the experimental knowledge 
of excited states. As there is a time difference of about 
10 years between the two compilations, we expect the 
TAT set to be more up-to-date in this respect. Above 
the highest know excited states, an integral over level 
densities has to be performed. The above mentioned 
differences in the level density treatment would also 
apply to the partition functions. 

4.2. A Global Comparison 

Following the discussion in §4.1.2 and 4.1.3, we 
expect a scattering of a factor 2-3 between the two 
rate sets due to the different treatments of gamma 
widths and level densities. On top of this, the sys- 
tematic deviations due to the different treatment of 
particle transmission coefficients, as discussed in sec- 
tion |4.1.l| , will enter. It is educational to perform 
an overall comparison. This is done in the following 
only for nuclei for which the ground state properties 
are identical in both sets, and for reactions that have 
identical Q-values. This limits the comparison to the 



vicinity of the valley of beta-stability and permits a 
more clearly defined comparison. It should, however, 
be mentioned that the usage of different mass mod- 
els for nuclei further from stability, entering reaction 
Q-values and ground state properties, introduces ad- 
ditional uncertainties. 

Hydrostatic oxygen burning is the first burning 
stage where a large number of reactions with unstable 
targets enter nucleosynthesis calculations. These are 
rates which are not often contained in the experimen- 
tal rate compilations. The temperature for this burn- 
ing stage is of the order 2 x 10 9 K. The next burning 
stage, silicon burning, approaches already a nuclear 
statistical equilibrium (NSE), where abundances are 
rather determined by chemical potentials (or bind- 
ing energies or reaction Q-values) and not by the in- 
dividual rates. Temperatures around 2 x 10 9 K are 
also of importance for explosive burning. At (ini- 
tially) higher temperatures in explosive environments, 
full NSE is usually attained and abundances depend 
again on masses or Q-values. In an expanding and 
cooling medium such equilibria freeze-out typically at 
temperatures of «2.5 x 10 9 K. It is here where indi- 
vidual reaction rates count and have the strongest 
influence. Thus, for various reasons, the tempera- 
ture range around 2.5 x 10 9 K is the most important 
for thermonuclear reaction rates of intermediate and 
heavy nuclei, and we chose this temperature for a gen- 
eral comparison. 

4-2.1. Comparison: Theoretical Rates 

In order to understand the differences, it is instruc- 
tive to take a look at global trends in the comparison 
of the TAT and WFHZ rates in Figures The ra- 
tios of the TAT/ WFHZ rates are plotted for a large 
number of nuclei for (n,7), (p,n), (a,n), and (a,p) re- 
actions. Let us start with the discussion of (n,7) rates. 
A temperature of 2.5xl0 9 K corresponds to thermal 
neutron energies of the order 0.25 MeV. For such en- 
ergies Figure m already indicates only small deviations 
of less than a factor of two for neutron transmission 
coefficients (somehwere between the dotted and short- 
dashed line) . This is smaller than the remaining scat- 
ter in the comparison of the (n,7) rates, due to the 
different treatment of other properties, like gamma 
transmission coefficients and nuclear level densities. 
Thus, we find a relatively flat deviation for the ratios 
of transmission coefficients with an overlying scatter 
of a factor 2-3, as expected from sections 1.1.2 and 
4.1.3 . This scatter is not specifically pronounced as a 
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function of nuclear structure and close to a random 
behavior. 

To avoid this general scatter, we will focus in the 
following on particle channels, as they result from the 
use of different optical potentials, where we expect 
more systematic differences. For this we choose as 
a next example in Figure |?j (p,n) reaction rates. A 
temperature of 2.5xl0 9 K corresponds in the Gamow 
window to an energy of «2.5 MeV protons. For such 
energies we see in Figure ^ strong deviations from 
unity and expect a maximum around Asa70-80. This 
is exactly what can be seen in Figure |^. We also see a 
minimum around j4w40-50 and again a rise to smaller 
mass numbers, as expected from Figure |J. 

Given the quite different values for, e.g., the proton 
transmission function (Fig. |^), it is perhaps surpris- 
ing that the rates (e.g., for (p,n) reactions, Fig. ^) dif- 
fer so little. This is because the (p,n) rates on heavy 
nuclei are only astrophysically important at temper- 
atures above about 3 x 10 9 K where the Gamow en- 
ergy is above 2 MeV. Figure f| is only for s- waves, but 
other partial waves do contribute. The difference in 
the proton transmission function ratio will decrease 
with increasing partial wave because the angular mo- 
mentum barrier will dominate and this barrier is the 
same in both approaches. At lower energies we ex- 
pect that the rates would differ more because of two 
reasons: The higher partial waves will be more sup- 
pressed and the difference in s-wave proton transmis- 
sion function ratio becomes slightly more pronounced. 
However, when going from 2 MeV to 1 MeV we would 
not expect much difference. Only at even lower en- 
ergies would an increase be noticeable but these are 
strongly hindered by the Coulomb barrier. 

In Figure || we display the rate comparison for 
(a,n) reactions. Again, as the neutron transmission 
coefficient ratio behaves already flat for these temper- 
atures, we expect essentially to see the behavior of 
alpha transmission coefficients. This is also exactly 
what can be observed (see the final paragraphs of 
|4.1.1| ). The ratio rises as a function of A, as also seen 
in Figure |5|, with the deviations from unity starting 
around A— 60. Figure ^ for the (a,p) rates compari- 
son combines the effects of Figure [7] with a maximum 
at A= 70-80 and Figure || with a continuous rise for 
A>60. We see deviations in this upper mass range 
which can amount to a factor 8. Up to A ss 55 we 
still see a smaller scatter of a factor 2-3. 

Although the above effects of utilizing different po- 
tentials can already be found for reactions with the 



target being in the ground state, they may have fur- 
ther consequences in the calculation of stellar rates in- 
volving thermally excited targets. As each thermally 
populated target state gives rise to similar deviations 
as described above, the temperature dependence of 
the rates will also be altered when comparing WFHZ 
and TAT rates. 

One additional difference should be mentioned here. 
TAT and WFHZ made different assumptions about 
isospin mixing, which is most important for Z=N nu- 
clei, and accounts for the well known inhibition of 
dipole transitions in self-conjugate (zero isospin) nu- 
clei. From the sparse experimental data available at 
the time (Toevs, 1971; Rogers et al. 1977; Dixon & 
Storey, 1977; Cooperman, Shapiro, & Winkler, 1977), 
Fowler and Woosley applied an empirical correction 
factor of 0.2 to the photon transmission function for 
these reactions. A less dramatic suppression of 0.5 
was assumed for (11,7) and (p,7) reactions into self- 
conjugate nuclei, and a factor of 0.67 for all photon 
transmission functions in \Z — N\ = 1 nuclei. There 
is less experimental justification for the latter. These 
correction factors were not included in the TAT rates 
nor in the TNH nucleosynthesis calculations. This 
accounts for the roughly 30% difference in the yield 
ratio for 44 Ti (§2). We defer discussion of these issues 
as they relate to 44 Ti synthesis until §5. 

4-2.2. Comparison: Theory vs. Experiment 

An important test of any theoretical result is how 
well it compares to experiment. The Hauser-Feshbach 
codes themselves are constructed with a wealth of 
experimental data (isotopic binding energies, ground 
state and excited state level energies and their as- 
sociated spins & parities, level densities, etc.). How 
this translates into reliable reaction rates can only 
be determined by direct comparison to published ex- 
perimental data. In this section we compare cross 
sections predicted by the two Hauser-Feshbach codes 
("CRSEC" for WFHZ & "SMOKER" for TAT) and 
reaction rates calculated from those cross sections, 
against their experimentally measured counterparts. 
The body of experimental reaction rate data was 
drawn from the literature up to the year 1992 (for 
a complete list of authors and reaction channels stud- 
ied, see Hoffman and Woosley (1992) or the www ad- 
dress previously listed). 

Wherever experimental rate information is avail- 
able, it should be used in preference to a theoreti- 
cal prediction. This was the principle motivation of 
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both the TAT and WFHZ groups in adopting the re- 
sults from the numerous reaction rate compilations of 
W. A. Fowler and his collaborators for reactions on 
light and intermediate mass nuclei (A §30). Beyond 
this the number of stable species (for each element), 
the large number of particle channels, and the impor- 
tant reactions that proceed through them, coupled 
with the experimental difficulties encountered when 
dealing with Coulomb barrier penetration on progres- 
sively heavier targets, conspire to limit the number 
of experiments that have been carried out on targets 
heavier than silicon. 

We begin by considering measured neutron capture 
cross sections. This will provide a valuable check on 
the reliability of the photon-transmission functions 
adopted by the TAT and WFHZ efforts. Shown in 
Figures 10 and |ll| are the ratios of the theoretical 
neutron capture cross sections (in mb and evaluated 
at 30 keV) divided by their recommended experimen- 
tal values (Bao & Kappeler 1987) vs. mass number 
A=Z+N. Where available, the error in the measured 
cross section divided by the cross section itself is in- 
dicated above and below unity (the dashed line) by 
error bars. Cross section ratios (WFHZ/EXP) are 
given by a filled diamond, (TAT/EXP) are given by 
an open square. Pairs of ratios reflecting reactions on 
a common target are coupled by a solid vertical line, 
the element of the target is listed above the largest 
ratio for a given comparison, with the mass num- 
ber given on the abscissa. Horizontal dotted lines 
are a factor of two above and below unity. Not 
shown in Figure [l^ are the cross section ratios for two 
targets, 26 Mg (WFHZ/EXP=16.7, TAT/EXP=17.6), 
and 31 P (WFHZ/EXP=4.5, TAT/EXP=4.7). For 
21 Ne the theoretical to experimental cross section ra- 
tio for WFHZ is 1.8, while for TAT it is 2.8 (off scale). 

We can calculate a statistical measure of the abil- 
ity of the Hauser-Feshbach codes to predict cross sec- 
tions (and reaction rates) over the range of available 
experimental data by determining the mean and stan- 
dard deviation of the value of the ratios for cross sec- 
tions predicted by theory divided by those measured 
by experiment. We restrict our comparisons to cross 
sections and reaction rates on targets with A > 28. 
In general the level density for reactions on lighter 
targets is too low for Hauser-Feshbach studies to be 
valid. Moreover, both groups have used experimental 
reaction rates for A < 28, we are testing our ability 
to model reactions on heavier targets. 

For the 30 keV cross sections (starting at 28 Si) 



shown in Figures [H] and [fl] the mean (and stan- 
dard deviation) of the cross section ratios predicted by 
SMOKER and CRSEC are 1.08 (0.64) and 0.91 (0.58) 
respectively. For the cross sections that are compared, 
both Hauser-Feshbach codes predict (n,7) cross sec- 
tions that are close to experiment nearly equally of- 
ten. On average SMOKER predicts a cross section 
that is larger than experiment, CRSEC predicts cross 
sections that are lower than experiment. The spread 
around the means is smaller for CRSEC, but not by 
much. The somewhat smaller deviations for CRSEC 
in this comparison with known experimental rates for 
stable nuclei also reflect to some extent the utiliza- 
tion of experimental level densities where available 
(see 4.1.3), i.e. predominantly for such nuclei where 
experimental cross sections exist. Therefore, in pre- 
dictions for unstable nuclei where such information is 
not available, one should probably expect (at best) 
accuracies similar to those resulting from SMOKER 
(where only theoretical predictions where used) . 

For charged-particle induced reactions on interme- 
diate mass nuclei, the experimental reaction rates im- 
plemented in the WFHZ reaction library are mostly 
drawn from the efforts of D.G. Sargood and his col- 
laborators (see §2 and Hoffman & Woosley 1992). 
We present comparisons for 63 experimentally mea- 
sured charged-particle induced reaction rates, sepa- 
rated into five reaction channels: (p,n), (p,7), (a,p), 
(a,n), and (a, 7). All comparisons will be made at 
one temperature, Tg = 3.0. 

Figures and 13 show ratios of (p,n) and (p,7) 



reaction rates that were calculated from the two the- 
oretical cross section codes (SMOKER and CRSEC) 
divided by experimentally measured reaction rates for 
protons incident on the targets listed by the element 
symbol in the figure and the mass number on the ab- 
scissa. In Figure [l^ the lightest reaction compared 
is for 27 Al(p,n) 27 Si. These figures are similar in all 
respects to Figures [lC] &: JTlj , but no error bars are 
given. Figures |lj and |l5[ show the same ratios for 
alpha- induced reactions. Statistical measures of the 
average error for Figures [w]- ^| are compiled in Table 
[|, which lists the number of theoretical reaction rates 
(or cross sections) compared to experiment, and the 
mean and standard deviations of these ratios for each 
reaction channel studied. 



For all reaction channels compared in Figures 10 
- each theoretical code predicted cross sections 
leading to reaction rates that agreed with experiment 
roughly an equal number of times, although for the 
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(p,n) rates the mean error for the TAT/EXP reaction 
rate ratio is closer to unity and the spread about the 
mean is less. For both codes, the error statistics are 
weighted heavily by the two (p,n) rates on 45 Sc and 
50 Ti, all other rates falling within the "factor of two" 
reliability lines. If these ratios are removed from the 
error analysis, the means (and standard deviations) 
for WFHZ/EXP and TAT/EXP reaction rate ratios 
then become 1.24 (0.32) and 0.85 (0.40) respectively 
with the mean again closer to unity for the more re- 
cently calculated rates, although the spread around 
each mean is now marginally better for the WFHZ 
set. This may reflect the inclusion of more and better 
known excited state data, and superior particle trans- 
mission functions inherent in the TAT treatment. 

For the (p,7) rates, the statistical errors indicate 
that CRSEC predicted reaction rates closer to unity, 
but again the results are weighted heavily by one 
problematic reaction rate. 39 K(p,7) 40 Ca has a WFHZ 
reaction rate to experiment ratio of 3.22, while the 
TAT ratio is much larger, 9.23 (off scale in Fig. 
fl3| ). Neither Hauser-Feshbach code predicted this 
rate with much accuracy, and if it is excluded from 
the error analysis the statistics improve considerably, 
with the mean (and standard deviation) for SMOKER 
and CRSEC being 1.33 (0.65) and 1.35 (0.81) respec- 
tively. These values are reported in Table [|. An- 
other reaction (not considered in these comparisons) 
is 41 K(p,a) 38 Ar, which is endoergic (has a positive 
"Q" value in the reverse reaction channel sense, i.e. 
decreasing charge or mass). This (p,a)-reaction rate 
(evaluated at Tg = 3.0), as predicted by both the 
SMOKER and CRSEC statistical codes, was under- 
estimated by a factor of three. Both of these cases are 
examples of possible deviations for reactions close to 
shell closures where level densities and gamma trans- 
mission coefficients undergo large changes. 

For the alpha-induced reaction channels, Figures 
pT| and |l5| show that nearly all of the comparisons 
fall within the "factor of two" reliability lines. The 
statistical errors for the (a,n) rate comparisons from 
each code are very nearly the same, while for the 
(a,p) comparisons they are only marginally differ- 
ent (Table |]). Finally, for the (0,7) channel, we 
compare only two experimentally measured reaction 
rates. For 34 S(a, 7) 38 Ar, the reaction rate ratios 
are WFHZ/EXP=1.35 and TAT/EXP=1.07. For 
42 Ca(a, 7) 46 Ti, the same numbers are 0.84 and 1.43 
respectively. Both codes predict these rates well 
within a factor of two uncertainty. 



It should be noted that for all reaction channels 
presented in Figures [lO] through [l5| and Table ||, the 
largest discrepancies are often confined to a few rates, 
particularly those either on (and especially leading 
into) targets with closed neutron shells. This is where 
the two different treatments of the nuclear level den- 
sity display the largest uncertainties. Nuclear shell 
corrections are very pronounced at or near magic 
nucleon numbers. Since reliable experimental data 
on radiation widths are lacking in these regions, it 
makes it difficult to disentangle the level density ef- 
fects from a possible problem in the prediction of the 
gamma transmission functions [although see also the 
discussion concerning the assumed nuclear potent ials 
(Equivalent Square Well vs. Woods-Saxon) in § 3.2 
and Figures HH . 

Reactions on targets in the vicinity of N=20 have 
received special attention by Sargood and his collabo- 
rators. Specific examples are the proton-induced reac- 
tions on 50 Ti, and 52 Cr. The 41 K(p,a) 38 Ar rate leads 
into this same closed shell. For the (a,p) channel reac- 
tions on both 48 Ti and 50 Ti, the statistical codes per- 
form poorly. Two exceptions are 34 S(a, 7) 38 Ar and 
54 Fe(a,n) 57 Ni, both were well calculated by the sta- 
tistical model codes. It should again be noted here 
that all of the experimental reaction rate information 
discussed above supplements the WFHZ reaction rate 
set which was utilized for the stellar yield predictions 
of §2, while TNH always made use of their theoretical 
TAT rates for intermediate mass nuclei. 

4.3. Future Improvements 

Based on the well-known statistical model code 
SMOKER (Thielemann, Arnould, & Truran 1987), 
an improved code for the prediction of astrophysi- 
cal cross sections and reaction rates in the statisti- 
cal model is presently developed. In the new code 
NON-SMOKER (Rauscher and Thielemann 1998), 
transmission coefficients for neutrons, protons and 7 
transitions are determined in the same way as previ- 
ously described in Section 4.1.1 . Optical potentials 
for a particles are treated in the folding approach of 
Satchler & Love (1979), with a parameterized mass- 
and energy-dependence of the real volume integral 
(Rauscher 1998). The mass- and energy-dependence 
of the imaginary potential is parameterized according 
to Brown and Rho (1981) and additionally includes 
microscopic and deformation information (Rauscher 
1998). 

The level density treatment has been recently im- 
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proved (Rauscher, Thielemann, & Kratz 1997). How- 
ever, the problem of the parity distribution at low 
energies remains. The new code includes this modi- 
fied description, but also accounts for non-evenly dis- 
tributed parities at low energies, based on the most re- 
cent findings within the framework of the shell model 
Monte Carlo method (Nakada & Alhassid 1997). 

Additionally, the included data set of experimen- 
tal level information (excitation energies, spins, par- 
ities) has been updated with the aid of the present 
version of the Evaluated Nuclear Data File (ENSDF; 
Firestone 1996), as well as the experimental nuclear 
masses (Audi and Wapstra 1995). For theoretical 
masses, currently the FRDM (Moller et al. 1995) is 
favored. Microscopic information needed for the cal- 
culation of level densities and a+nucleus potentials 
are also taken from the FRDM, as well as experi- 
mentally unknown ground state spins (Moller, Nix, 
& Kratz 1997). 

Finally, isobaric analog states T y = T K + 1 = 
T g s - + 1 are explicitly considere d in the new code 
( Rauscher and Thielemann 1996 ) . This replaces the 
empirical constant suppression factor of the 7-width 
used in the WFHZ rates for self-conjugate nuclei. 

For comparison, we give in Table || a select set 
of preliminary reaction rates obtained with the code 
NON-SMOKER (Rauscher and Thielemann 1998), 
along with their counterparts in the TAT and WFHZ 
rate libraries. These rates are identified (in the next 
section) as having influenced the observed discrepan- 
cies in the stellar yields for the hybrid models of § 



5. Nuclear Origin of the Nucleosynthetic Dif- 
ferences 

The nuclei which show greater than 20% deviations 
in cither the 15 M Q or the 25 M Q star are: 33 S, 40 Ar, 
40 K, 44 ' 46 Ca, 45 Sc, 50 V, and 70 Zn. Each of these will 
now be discussed. Quoted ratios of stellar yields and 
reaction rates here will always be given as those of 
TAT divided by WFHZ. The related reaction rates 
are given m Tables | and p. 

The nucleus 33 S is predominantly made in explo- 
sive oxygen burning (Tg ~ 3.5) and explosive neon 
burning (T 9 ~ 2.5) (Thielemann & Arnett 1985, 
WW95). The main production and destruction re- 
actions are 32 S(n,7) 33 S and 33 S(p,7) 34 Cl respectively. 
30 Si(a,n) 33 S (and its inverse) can also compete, de- 
pending on the neutron and a-particle abundances 



and temperature. It can even be the dominant chan- 
nel for production or destruction. In the 15 M Q and 
25 M stars, the 33 S yield ratio (Y(TAT) /Y(WFHZ), 
Figures |l]and||) differed by factors of 1.23 and 1.33, 
respectively. At a temperature of Tg = 3.5 as shown 
in Table |5|, the ratio of the production rates for 
32 S(n,7) 33 S is up by a factor of 1.3 and the destruction 
rate 33 S(p,7) 34 Cl is up by a factor of 2.0. However, 
33 S(n,a) 30 Si is down by 0.7. This combination, in 
addition to the particle fluxes at these high temper- 
atures, leads to the observed enhancement. It seems 
the (a,n) reaction works more in favor of destruction. 

40 Ar and 40 K are the result of neutron capture dur- 
ing hydrostatic helium, carbon, and to a lesser extent, 
neon burning. Since these two nuclei are linked by a 
common reaction flow, the reasons for their observed 
discrepancies in the stellar yields are likewise coupled. 
The Ar yield ratio is down by factors of 0.73 and 
0.62 while 40 K is up by a factor of 1.21 and 1.33, 
respectively, in the 15 and 25 M Q stars. 

In helium and carbon burning the neutron sources 
are 22 Ne(a,n) 25 Mg and 13 C(a,n) 16 0. Which domi- 
nates depends on the stellar conditions and current 
state of evolution. At carbon burning temperatures 
(~ 1 x 10 9 K), 40 K is produced by 39 K(n, 7 ) 40 K and de- 
stroyed by 40 K(n,a) 37 Cl or 40 K(n,p) 40 Ar. While the 
ratio of the production reaction rates is slightly en- 
hanced (by a factor 1.2), the two destruction reaction 
ratios are suppressed by factors of 0.77 and 0.4. This 
leads to the enhancement of 40 K by factors of 1.21 and 
1.33 seen in the 15 and 25 Mq stars respectively. Sim- 
ilarly the ratio of the production rate, 39 Ar(n,7) 40 Ar, 
is near unity, but the production is reduced by a fac- 
tor 0.4 in the second channel 40 K(n,p) 40 Ar mentioned 
above. This leads to the reduced production by fac- 
tors of 0.73 and 0.62 seen in Figures [l] and ^ . 

Ca is a product of neutron capture in hydrostatic 
helium, carbon and neon burning. In the 15 and 25 
M models the yield ratios are 1.8 and 1.88 respec- 
tively. In the chain 45 Ca(n,7) 46 Ca(n,7) 47 Ca the ratio 
of the producing (n,7) reaction rates is higher by a 
factor of 1.6, while the destruction reaction rate ratio 
is near unity. This accounts for most of the alteration 
in the yield. 

A special note should be made here in reference to 
the implementation of experimental (n,7) cross sec- 
tion data that lead to three of the reaction rates given 
in Table ||. This accounts for some of the differences 
seen in the neutron capture rates on 32 S, 39 K, and 
46 Ca which affect the production of 33 S, 40 K and 46 Ca 
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respectively. In the WFHZ rate set, the leading pa- 
rameter in the fit to the theoretical (11,7) rate ("A" in 
eq. 3, WFHZ, 1978) has been modified to give a re- 
action rate at 30 keV that agrees with the suggested 
30 keV cross section of Bao & Kappeler (1987), while 
keeping the temperature dependence of the original 
theoretical (n,7) rate. The TAT set also made use of 
these experimental rates, assuming however a domi- 
nant s-wave behavior which leads to a constant rate as 
a function of temperature. Therefore both implemen- 
tations only agree at about 30 keV (= 3 x 10 8 K), with 
the result that the TAT rate under(over)estimates the 
(n,7) rate at lower(higher) temperatures, consistent 
with the 1/v behavior of neutron capture rates, al- 
though the deviation in most cases is not large, often 
less than a factor of three, over a wide range in tem- 
perature. This leads to the following different neu- 
tron capture rates at the relevant temperatures in the 
WFHZ and TAT rate sets, respectively: For 32 S at 
T 9 =3.5 the (n, 7 ) rates are 5.4 x 10 5 and 6.9 x 10 5 , for 
39 K at T 9 =1.0 they are 1.4x 10 6 and 1.7x 10 6 , and for 
46 Ca at r 9 =1.0 7.8 x 10 5 and 7.9 x 10 5 (all rates given 
in cm 3 s -1 mole -1 ). Thus, the overproduction of 33 S 
and 40 K, respectively, in the TAT models as compared 
to the WFHZ models would be less pronounced with 
the same implementation of the experimental rates. 

44 Ca is dominantly produced as 44 Ti in the a-rich 
freeze-out phase of explosive silicon burning. It is 
thus the product of the a-capture chain along self- 
conjugate a-nuclei 28 Si, 32 S, 36 Ar, 40 Ca and 44 Ti. 
In Table |^ we give the stellar (a, 7) reaction rates 
on self-conjugate nuclei from Mg to Ti (all evaluated 
at Tg = 2.5) as derived from the statistical model 
codes CRSEC, SMOKER, and NON-SMOKER, re- 
spectively. 

For the (a, 7) rates up to and including 24 Mg(a, 7) 28 Si, 
the same experimental information (CF88) has been 
used in both rate sets. The reason for the larger 
SMOKER rates is the neglect of isospin suppression 
effects for self-conjugate nuclei that were included in 
the CRSEC rates. This effect is apparent in the 
new NON-SMOKER rates, which, at least up to 
36 Ar, are in good agreement with those of CRSEC. 
The 24 Mg(a, 7) 28 Si reaction cross section has been 
measured (CF88) and the stellar rate (3.7 cm 3 s _1 
mole -1 at T g =2.5) is in very good agreement with 
the new NON-SMOKER result (2.8 cm 3 s" 1 mole- 1 ). 
In the case of 40 Ca(a, 7) 44 Ti, experimental data is 
also available (Cooperman, Shapiro & Winkler 1977) 
but was not used in the calculation of this cross sec- 



tion by either TAT or WFHZ. The experiment sug- 
gests a value in the temperature range of interest 
(Tg w 2.5) which is about 4 times lower than the 
CRSEC rate, in excellent agreement with the NON- 
SMOKER prediction, and in considerable disagree- 
ment with SMOKER (Wiescher, private communica- 
tion). All of this is encouraging, and indicates an 
area where modern laboratory measurements might 
lead to a more reliable determination of some very 
important reaction rates. 

It should be noted for the reactions that directly 
impact 44 Ti synthesis that their exact values will 
not affect abundances as long as temperatures are 
high enough to sustain an (a, 7) — (7, a)-equilibrium, 
wherein the abundances are only determined by Q- 
values and partition functions, but the different rates 
will partially enter during the charged particle freeze- 
out around 3xl0 9 K, resulting in the enhancement of 
the 44 Ti yield ratio by factors of 1.28 and 1.31, respec- 
tively, for the 15 and 25 M models. This is particu- 
larly apparent on consideration of the large disparity 
between the TAT and WFHZ rates for the dominant 
production reaction ( 40 Ca(a, 7) 44 Ti, which differs by 
a factor of nearly 6) , and the dominant destruction re- 
action [ 44 Ti(a,p) 47 V, which is identical between the 
two sets and also shows no change in the preliminary 
NON-SMOKER calculation, see Tables | and § . This 
is a nice example of where large reaction rate differ- 
ences cause much smaller abundance deviations due 
to equilibria at high-temperatures in explosive burn- 
ing and only enter during the final freeze-out phase. 

Several other reaction rates of importance to 44 Ti 
synthesis have been identified by The et al. (1998), 
who used the TAT rate set in their studies. One such 
reaction not previously mentioned is 45 V(p,7) 46 Cr 
which is on the border of a QSE cluster that af- 
fects 44 Ti. At temperatures where the material is 
falling out of equilibrium (Tg w 2.5) this rate be- 
comes important, and the TAT rate is in severe dis- 
agreement with the WFHZ rate, with a reaction ratio 
of TAT/WFHZ=0.016. This had no impact on the 
present calculations, because Cr was not included in 
our networks. The et al. also point out that for a com- 
position with a non-zero neutron excess (rj > 0.0004), 
this reaction rate does not affect 44 Ti synthesis. 

45 Sc is made in two complementary ways via the 
radioactive progenitor 45 Ti in explosive oxygen and 
silicon burning, and directly during neon and car- 
bon burning. In the 25M star, and using the TAT 
rate set, 45 Sc was made in equal proportions as itself 
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and radioactive 45 Ca. When using the WFHZ rare 
set, it was made (predominantly) as itself. In the 
15 M-q star the production via 45 Ti in explosive oxy- 
gen and silicon burning dominates. The production 
is reduced by a factor of 0.85 in the yield comparison 
TAT/WFHZ. This occurs via the reaction sequence 
42 Ca(a,n) 45 Ti(n,p) 45 Sc and 43 Ca(p, 7 ) 44 Sc(p,7) 45 Ti. 
In the first sequence, the (n,p)-reaction is enhanced 
by a factor 1.4 (i.e. enhanced 45 Ti destruction). The 
first rate in the proton capture sequence, which con- 
trols the flow, is reduced by a factor 0.9. In combina- 
tion this lead to the reduction of a factor of 0.85. 

50 V is made in explosive and hydrostatic neon 
burning, coupled with some neutron processing. The 
destruction rate 50 V(n,7) 51 V is reduced by a factor of 
0.75, thus leading to an enhancement of factors 1.25 
and 1.28 respectively for the 15 and 25 M stars. 

The heavy nucleus 70 Zn is made by neutron cap- 
ture in He, C, and Ne-burning, i.e., the s-process. It's 
abundance variation can be explained by changes in 
neutron capture cross sections. The reaction chain 
69 Zn(n,7) 70 Zn(n,7) 71 Zn has the first rate increased 
by a factor of 1.97, the second by 1.7. This lead to a 
net gain of a factor 1.22. 

6. Summary and Conclusions 

We have reviewed the theoretical nuclear reaction 
rates currently used for modern calculations of stellar 
nucleosynthesis in the intermediate mass range (mag- 
nesium through krypton) and have explored the sen- 
sitivity of standard calculations to those rates. Both 
sets, the one of Thielemann et al. (TAT) and the 
one of Woosley et al. (WFHZ), are derived from 
Hauser-Feshbach theory. We discussed the underly- 
ing assumptions in each rate set, providing much pre- 
viously unpublished information. We compared them 
to one another and to experiment, elucidating the dif- 
ferences, and explored the nucleosynthesis obtained 
using each rate set in otherwise identical models for 
the evolution and explosion of 15 and 25 M super- 
novae. 

Overall, the two rate sets are similar. Typical dif- 
ferences at astrophysically interesting temperatures 
are less than a factor of two (Figures |-|). There 
are individual cases, however, where the difference 
exceeds a factor of 10. Some of the larger differ- 
ences occur for reactions where scarce experimen- 
tal information is available and different assumptions 
were made regarding the photon transmission func- 



tion, for example, (a, 7) reactions on Z = N nuclei. 
Different assumptions were also made about the par- 
ticle transmission functions, nuclear partition func- 
tions, and level densities. More modern and complete 
data used in the TAT rates makes them superior in 
cases where the partition function is important. As- 
sumptions regarding the particle transmission func- 
tions seem less important (see, however, Rauscher 
1998 and Somorjai et al. 1998 with respect to al- 
pha transmission functions for heavy nuclei). WFHZ 
used an equivalent square well with empirical reflec- 
tion factors; TAT used a more detailed optical model. 
Given the quite different values for, e.g., the neutron 
and proton transmission function (Figures ^ and |j), 
it is perhaps surprising that the rates (e.g., for (p,n) 
reactions, Fig. ?]) differ so little. This is because 
the relevant temperatures for explosive burning are 
high. For incident particles in the Gamow window, 
the deviations in the particle transmission functions 
are typically smaller than a factor of two. In addition, 
higher partial waves (not shown in Fig. contribute. 
A comparison of rates at a lower temperature would 
have revealed larger discrepancies. 

Compared to experiment, both sets of theoretical 
rates give similar agreement, typically to a factor of 
two (Figures |l0| - |l^) . The standard deviations be- 
tween the two theoretical sets and data for (n,7) cross 
sections - which measure how well the photon trans- 
mission function is calculated - are virtually identical, 
0.58 for WFHZ/EXP and 0.64 for TAT /EXP (§4.2.2) 
and both sets predict a mean that is equally close to 
experiment, TAT/EXP is higher (1.08), WFHZ/EXP 
is lower (0.91). For charged particles - which are more 
sensitive to particle transmission functions - the TAT 
calculations are slightly superior, with standard devi- 
ations over the entire range of experimentally know 
reaction rate data being 0.81 and 0.65 [for (p,7) rates] 
and 0.62 and 0.47 [(p,n) rates] for WFHZ/EXP and 
TAT/EXP respectively. Other comparisons are given 
in Table | and §4.2.2. 

In summary, the two rate sets have comparable 
merit when compared to experiment. All the au- 
thors of this paper agree that the new rate set - the 
"NON-SMOKER" set, will be preferable to both TAT 
and WFHZ and, when they become available, will be 
adopted by both groups for future work. 

When the two current rate sets are included in oth- 
erwise identical stellar models we find that the nu- 
cleosynthesis, again with rare and occasionally inter- 
esting exceptions, is not greatly changed. Individual 
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exceptions are discussed in §5.1. For example, only 
about a dozen (out of 70) stable isotopes in the mass 
range 12 to 70 have nucleosynthesis that differ by over 
20% in two supernovae of 15 M that use the same 
rate for 12 C(a,7) 16 (Table |J). It can, however, be 
noticed that most of these isotopes - with one excep- 
tion 44 Ti - are products of hydrostatic burning where 
individual reaction flows are governed by the cross 
sections involved. Nevertheless, none differ by more 
than a factor of 1.7. Given the significantly larger 
differences that exist in individual reaction rates, one 
may wonder at the robust nature of the final nucle- 
osynthesis. We see three major causes. 

First, as the star burns and becomes hotter, the 
nuclear flow follows the valley of beta-stability mak- 
ing heavier nuclei as it goes. In doing so, it follows the 
path of least resistance - those reactions having the 
largest cross section for a nucleon or a-particle react- 
ing with a given nucleus. These large cross sections 
are reasonably well replicated by any calculation, nor- 
malized to experiment, that treats the Coulomb bar- 
rier and photon transmission function approximately 
correctly. Large differences may exist in rate fac- 
tors for reactions that are in competition, especially a 
small channel in the presence of a large one, but these 
small channels are frequently negligible, at least for 
the major abundances while they can cause larger dif- 
ferences when one is interested in the abundances of 
trace isotopes. 

If one is however interested in accurate abundance 
predictions resulting from these smaller flows in hy- 
drostatic burning stages, this can in most cases only 
be obtained by improving the reliability of the cross 
sections (and reaction rates) that determine these 
weak flows on light and intermediate mass nuclei. As 
new experimental information becomes available, a 
continuous improvement is therefore highly warranted 
(see e.g. the most recent compilation of experimental 
rates, NACRE, Angulo et al. 1999). 

Second, beyond oxygen burning, which is to say for 
nuclei heavier than calcium, nucleosynthesis increas- 
ingly occurs in a state of full or partial nuclear statis- 
tical equilibrium. There the abundances are given by 
binding energies and partition functions. So long as 
the "freeze-out" is sufficiently rapid, individual rates 
are not so important. 

Third, the reaction rates varied here were only 
those theoretical values from Hauser-Feshbach calcu- 
lations for intermediate mass nuclei, i.e., nuclei heav- 
ier than magnesium. The really critical reaction rates 



are, for the most part, those below magnesium. These 
reactions, like e.g., 12 C(a,7) 16 0, govern the energy 
generation, major nucleosynthesis, and neutron expo- 
sure in the star. The rest are perturbations on these 
dominant flows. 

This is not to say, however, that the nuclear and 
stellar details of heavy element synthesis are now well 
understood. Differences in the stellar model may ac- 
count, not just for 20% variation, but orders of mag- 
nitude (Table ^; column TNH). That is, uncertainty 
in stellar physics - especially the treatment of con- 
vection and how it is coupled (or not coupled) to the 
nuclear network - accounts for most of the differences 
in current nucleosynthesis calculations - provided such 
calculations use the same nuclear reaction rates below 
magnesium. 

Even in a perfect stellar model though, there will 
still be interesting nuclear physics issues. Stellar nu- 
cleosynthesis is becoming a mature field rich with di- 
verse and highly detailed observational data. The 
"factor of two" accuracy that was adequate in the 
past may not do justice to the observations of the fu- 
ture. There are many individual cases where the nu- 
clear physics uncertainty is still unacceptably large. 
We would close by just pointing out two of them. 

The suppression of radiative capture reactions into 
self-conjugate (isospin zero) nuclei is very uncertain. 
Past Hauser-Feshbach calculations have adopted em- 
pirical factors for this suppression (§4.2.1). Alpha- 
capture reactions, like 24 Mg(a, 7) 28 Si, 28 Si(a, 7) 32 S, 
44 Ti(a, 7) 48 Cr, are very important to nucleosyn- 
thesis in oxygen and silicon burning. The reaction 
40 Ca(a, 7) 44 Ti also directly affects the synthesis of 
44 Ti. Modern accurate determinations of most of the 
reaction rates are missing (as well as (p,7) reactions 
into the same nuclei). Measurements here would be 
most welcome. 

The Hauser-Feshbach rates are also only as good 
as the local experimental rates to which the necessary 
parameters of the calculation are calibrated. In that 
regard we would point out the near absence of charged 
particle reaction rate data for A z 70. Charged par- 
ticle reactions are important, especially on unstable 
nuclei, at significantly higher atomic weights in the re- 
process (Hoffman, Woosley & Qian 1997) and in the 
p-process (Rayet et al. 1995). 
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8 16 24 32 40 48 56 64 72 

number of nucleons 

Fig. 1. — Comparison of nucleosynthesis using TAT 
rates and WFHZ rates (for A > 28]. Plotted is the 
ratio of the final stellar yield (M ) for a 15 M super- 
nova evolved from the main sequence and calculated 
using TAT rates divided by the yield calculated in an 
identically derived star using WFHZ rates. Below A 
= 28 the rate sets were identical. All models used a 
rate for 12 C(a, 7 ) 16 from Caughlan & Fowler (1988) 
multiplied by 1.7. 
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Fig. 2. — The same as Fig. |l| but for a 25 M super- 
nova. 
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Fig. 3. — Ratios of neutron transmission functions 
obtained with the TAT potentials and the equivalent 
square well description of WFHZ for a set of stable 
nuclei with varying mass number A. Shown are the 
ratios for s-wave neutrons. Different line styles denote 
different bombarding energies of the target-projectile 
system in the center of mass frame: 0.01 MeV (solid), 
0.1 MeV (dotted), 1.0 MeV (short dashes), 2.0 MeV 
(long dashes), 5.0 MeV (dot - short dash), 10.0 MeV 
(dot - long dash), 15.0 MeV (short dash - long dash), 
20.0 MeV (solid). 
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Fig. 4. — Same as Figure || but for ratios of s-wave 
proton transmission functions. The line styles denote 
different energies: 1.0 MeV (solid), 2.0 MeV (dotted), 
5.0 MeV (short dashes), 10.0 MeV (long dashes), 15.0 
MeV (dot - short dash), 20.0 MeV (dot - long dash). 
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Fig. 5. — Same as Figure [| but for ratios of s-wave 
alpha-particle transmission functions. Line styles de- 
note different energies: 5.0 MeV (solid), 10.0 MeV 
(dotted), 15.0 MeV (short dashes), 20.0 MeV (long 
dashes) . 
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Fig. 6.— Ratios of TAT over WFHZ (11,7) rates 
evaluated at T 9 = 2.5. 
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Fig. 7.— Ratios of TAT over WFHZ (p,n) rates eval- 
uated at T 9 = 2. 5. The extreme values at A=64 and 
70 (0.1 and 20 respectively), give an indication of the 
differences that arise when calculating theoretical re- 
action rates using nuclear data far from stability. The 
two cases are for (p,n) rate ratios on targets 64 Ga and 
70 Se (each 5 and 4 units from stability towards the 
proton drip line) which have (p,n) Q-values that dif- 
fer by 97 and 540 keV (respectively) in the two rate 
sets. Also, in the first case the spin of the compound 
nucleus differs by 1, in the second case it differs by 
1 for the compound nucleus and 2 for the final nu- 
cleus. This leads to different channel openings due to 
Q-values and angular momentum barriers and shows 
how crucial mass and ground state properties are far 
from stability. 
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8.— Ratios of TAT over WFHZ (a,n) rates eval- 
3d at T 9 = 2.5. 



29 



10 



.1 — 



20 30 40 



50 60 
A 



70 80 90 



Fig. 9.— Ratios of TAT over WFHZ (a,p) rates eval- 
uated at T g = 2.5. 
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Fig. 10. — 30 keV neutron-capture cross section ra- 
tios vs. mass number (20< A < 50); Theory (WFHZ: 
filled diamonds, TAT: open squares) divided by ex- 
periment. Error bar values (error/cross section) are 
shown for each comparison ratio. Ratios for (11,7) 
cross sections on common targets are connected by a 
solid vertical line. For reactions proceeding on tar- 
gets with a common mass number, the second pair is 
circled (eg. 40 Ca). 
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Fig. 12. — Proton-induced (p,n) reaction rate ratios 
vs. mass number A; Theory (WFHZ: filled diamonds, 
TAT: open squares) divided by experiment. Ratios 
for (p,n)-reactions on common targets are connected 
by a solid vertical line. For reactions proceeding on 
targets with a common mass number, the second pair 
is circled (eg. 48 Ti). 
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Fig. 13. — Same as Figure [l^ but for (p,7) reaction 
ratios. 
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Fig. 14. — Same as Figure [12] but for (a,p) reaction 
ratios. 
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Fig. 15. — Same as Figure [12] but for (a,n) reaction 
ratios. 
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Nuclear Reaction Networks Employed 
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15 M Stellar Yields (in M ) 12 C(a, 7)=2.35*CF88 
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1.00 


6.71(-16) 




nse-Ia-MCh 


45 Sc 


3.55(-6) 


2.91(-6) 


1.22 


4.66(-8) 


0.016 


a,xO,xSi,C,Ne,i/ 


46 Xi 


1.81(-4) 


1.90(-4) 


0.95 


2.72(-6) 


0.014 


xO,Ia-det 


47 Xi 


1.02(-5) 


1.08(-5) 


0.94 


4.67(-6) 


0.432 


xO,xSi,Ia-det 


48 Xi 


1.56(-4) 


1.72(-4) 


0.91 


1.27(-4) 


0.738 


xSi,Ia-det 


49 Ti 


1.49(-5) 


1.18(-5) 


1.26 


4.15(-6) 


0.352 


xSi 


50 Ti 


2.77(-6) 


2.83(-6) 


0.98 


4.59(-10) 




nse-Ia-MCh 


50y 


1.06(-7) 


1.34(-7) 


0.79 


3.51(-10) 


0.003 


C,Ne,xNc,xO 


Sly 


4.36(-5) 


5.90(-5) 


0.74 


1.01(-5) 


0.171 


a,Ia-det,xSi,xO,z/ 



38 



Table 2 — Continued 



A z 


WFHZ 


TAT 


WFHZ /TAT 


TNH 


TNH/TAT 


process(es) 


50 Cr 


2.31( 


-4) 


2.83( 


-4) 


0.81 


3.85(-5) 


0.136 


xSi,xO,a,Ia-det 


52 Cr 


1.37( 


-3) 


1.43( 


-3) 


0.96 


8.24(-4) 


0.576 


xSi,a,Ia-det 


53 Cr 


1.70( 


-4) 


1.83( 


-4) 


0.93 


8.92(-5) 


0.487 


xO,xSi 


54 Cr 


8.08( 


-6) 


8.06( 


-6) 


1.00 


l.ll(-lO) 


— 


nse-Ia-MCh 


55 Mn 


1.12( 


-3) 


1.11( 


-3) 


1.00 


3.39(-4) 


0.305 


Ia,xSi,f 


54 Fe 


6.94( 


-3) 


7.31( 


-3) 


0.95 


3.51(-3) 


0.480 


Ia,xSi 


56 Fe 


9.65( 


-2) 


9.61( 


-2) 


1.00 


1.30(-1) 


1.353 


xSi,Ia 


57 Fe 


3.45( 


-3) 


3.44( 


-3) 


1.00 


4.64(-3) 


1.349 


xSi,Ia 


58 Fe 


1.58( 


-4) 


1.60( 


-4) 


0.99 


3.82(-10) 




He(s),nse-Ia-MCh 


59 Co 


3.00( 


-4) 


2.81( 


-4) 


1.07 


1.36(-4) 


0.484 


He(s),a,Ia,f 


58 Ni 


4.22( 


-3) 


4.34( 


-3) 


0.97 


6.64(-3) 


1.529 


a 




2.50( 


-3) 


2.67( 


-3) 


0.94 


3.13(-3) 


1.172 


a,He(s) 


61 Ni 


1.52( 


-4) 


1.49( 


-4) 


1.02 


1.46(-4) 


0.980 


a,Ia-det,He(s) 


62 Ni 


7.11( 


-4) 


8.91( 


-4) 


0.80 


1.00(-3) 


1.122 


a,He(s) 


64 Ni 


4.36( 


-5) 


4.59( 


-5) 


0.95 


9.17(-16) 




He(s) 


63 Cu 


4.01( 


-5) 


3.16( 


-5) 


1.27 


9.56(-7) 


0.030 


C,Ne 


65 Cu 


1.27( 


-5) 


1.42( 


-5) 


0.89 


7.69(-7) 


0.054 


He(s) 


64 Zn 


2.84( 


-5) 


2.67( 


-5) 


1.06 


1.41(-5) 


0.528 


j/-wind,a;,He(s) 


66 Zn 


5.56( 


-5) 


5.78( 


-5) 


0.96 


1.47(-5) 


0.254 


He(s), a, nse-Ia-MCh 


67 Zn 


3.27( 


-6) 


3.54( 


-6) 


0.92 


1.94(-8) 


0.005 


He(s) 


68 Zn 


1.63( 


-5) 


1.75( 


-5) 


0.93 


6.35(-9) 




He(s) 


70 Zn 


2.81( 


-7) 


3.12( 


-7) 


0.90 


3.19(-21) 




He(s) 


69 Ga 


2.35( 


-6) 


2.23( 


-6) 


1.05 


6.47(-13) 




He(s) 


71 Ga 


1.21( 


-6) 


1.30( 


-6) 


0.93 


2.29(-20) 




He(s) 


70 Ge 


7.32( 


-6) 


5.53( 


-6) 


1.32 


5.13(-15) 




He(s) 



a Nucleosynthetic process(es) responsible for the production where "BB" stands for the big 
bang (e.g. Walker et al. 1991); "L*" stands for low-mass stars (M< 8 M e.g. Renzini & 
Voli 1981) or in some cases AGB stars (Sackmann & Boothroyd 1992); 'V is the neutrino 
process (e.g. Woosley et al. 1990); "f-wind is the neutrino-driven wind from young neutron 
stars in Type II and lb supernovae (e.g. Hoffman, Woosley, & Qian 1997); "a" is the a- 
rich freeze out from nuclear statistical equilibrium (e.g. Woosley, Arnett, & Clayton 1973); 
"Novae" indicates classical novae (e.g. Woosley 1986); "la" is Type la supernovae (Nomoto, 
Thielemann, & Yokoi 1984 for ordinary la's which are assumed to ignite at the Chandrasekhar 
mass; Woosley & Eastman 1995 for the inner neutron-rich regions of ordinary la's- here called 
"nse-Ia-MCh" ; Woosley & Weaver 1994 for sub-Chandrasekhar detonation models "Ia-He- 
det"). All other entries refer to burning processes in massive stars, H, He, C, Ne, O, Si- 
responsible for producing the isotope, with a prefix "x" indicating explosive nucleosynthesis. 
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Table 3 

Radioactive Stellar Yields (in M ) in the 15 M Supernovae 



A Z 


WFHZ 


TAT 


WFHZ/TAT 


TNH 


TNH/TAT 


process (es) a 


22 Na 


8.10(-8) 


9.68(-8) 


0.84 


3.98(-8) 


0.411 


Novae 


26 Al 


3.07(-5) 


2.92(-5) 


1.05 


2.68(-6) 


0.092 


j/,xNe,H 


44 Ti 


6.06(-5) 


7.40(-5) 


0.82 


7.19(-5) 


0.972 


a 


56 Ni 


8.11(-2) 


8.06(-2) 


1.01 


1.30(-1) 


1.613 


a,nse-Ia-MCh 


60p e 


3.35(-6) 


5.51(-6) 


0.61 


1.25(-17) 




Ne,xNe,xHe 


60 Co 


5.63(-6) 


2. 60 (-6) 


2.17 


4.14(-14) 




Ne,xNe,xHe 



a see footnote table § 
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Table 4 

Statistical Errors of Reaction Ratios 





# 




&WFHZ 


XT AT 


OTAT 


( n ,7) 


62 


0.91 


0.58 


1.08 


0.64 


(P, n ) 


12 


1.47 


0.62 


0.95 


0.47 


(P>7) 


20 


1.35 


0.81 


1.33 


0.65 


(«,P) 


12 


0.81 


0.35 


0.73 


0.45 


(a,n) 


16 


0.85 


0.39 


0.86 


0.23 


(a,7) 


2 


1.10 


0.36 


1.25 


0.25 



a x and a indicate the arithmetic mean and stan- 
dard deviation for the stated number (#) of cross 
section [(n,7)] and reaction rate (all other channel) 
ratios (theory/experiment). 
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Table 5 

Comparison of stellar reaction rates in cm 3 s _1 mole -1 from the WFHZ, TAT, and 

NON-SMOKER reaction rate libraries. 



Reaction 


T 9 


WFHZ 


TAT 


TAT/ WFHZ 


NON-SMOKER 


32 S(n, 7 ) 33 S a 


3.5 


5.4(5) 


6.9(5) 


1.28 


8.4(5) 


33 S(p, 7 ) 34 Cl 


3.5 


5.6(3) 


1.1(4) 


1.96 


6.2(3) 


33 S(n,a) 30 Si 


3.5 


5.1(7) 


3.4(7) 


0.67 


3.4(7) 


39 K(n, 7 ) 40 K a 


1.0 


1.4(6) 


1.7(6) 


1.21 


1.8(6) 


40 K(n,a) 37 Cl 


1.0 


1.1(7) 


8.5(6) 


0.77 


7.6(6) 


40 K(n,p) 40 Ar 


1.0 


4.5(6) 


1.8(6) 


0.40 


1.8(6) 


45 Ca(n, 7 ) 46 Ca 


1.0 


1.6(6) 


2.5(6) 


1.56 


1.7(6) 


46 Ca(n, 7 ) 47 Ca a 


1.0 


7.8(5) 


7.9(5) 


1.01 


4.7(5) 


44 Ti(a,p) 47 V 


2.5 


9.2(-3) 


9.6(-3) 


1.04 


9.6(-3) 


45 Ti(n,a) 42 Ca 


3.5 


6.4(6) 


9.9(6) 


1.55 


1.0(7) 


45 Ti(n,p) 45 Sc 


3.5 


1.7(8) 


2.4(8) 


1.41 


2.5(8) 


43 Ca(p, 7 ) 44 Sc 


3.5 


3.9(4) 


3.5(4) 


0.90 


3.2(4) 


44 Sc(p, 7 ) 45 Ti 


3.5 


1.2(4) 


2.5(4) 


2.08 


1.6(4) 


50 V(n, 7 ) 51 V 


3.5 


2.8(6) 


2.1(6) 


0.75 


3.9(6) 


69 Zn(n, 7 ) 70 Zn 


1.0 


6.6(6) 


1.3(7) 


1.97 


1.1(7) 


70 Zn(n, 7 ) 71 Zn 


1.0 


1.4(6) 


2.4(6) 


1.71 


2.6(6) 



a Actual rate used by TAT and WFHZ is based on Bao & Kappeler 1987 
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Table 6 

Comparison of stellar (0,7) reaction rates (in cm 3 s _1 mole -1 ) on self-conjugate nuclei 

DERIVED FROM THE CRSEC, SMOKER, AND NON-SMOKER STATISTICAL model CODES. 



Reaction 


T 9 


CRSEC 


SMOKER 


SM/CR 


NON-SMOKER 


24 Mg(a, 7 ) 28 Si 


2.5 


1.1(1) 


3.7(1) 


3.36 


2.8(0) 


28 Si(a, 7 ) 32 S 


2.5 


3.8(-l) 


4.5(0) 


11.71 


4.6(-l) 


32 S(a, 7 ) 36 Ar 


2.5 


2.3(-l) 


8.9(-l) 


3.87 


9.4(-2) 


36 Ar(a, 7 ) 40 Ca 


2.5 


3.4(-3) 


2.4(-2) 


7.06 


1.6(-3) 


40 Ca(a, 7 ) 44 Ti 


2.5 


1.5(-2) 


7.8(-2) 


5.70 


3.3(-3) 


44 Ti(a, 7 ) 48 Cr 


2.5 


3.3(-6) 


1.6(-4) 


48.48 


1.8(-5) 
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